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The two-dimensional repulsive Hubbard model has been investigated by a variety of methods, 
from small to large U. Superconductivity with d-wave symmetry is consistently found close to half 
filling. After a brief review of the various methods a variational many-electron state is discussed 
in more detail. This trial state is a natural extension of the Gutzwiller ansatz and provides a 
substantial improvement thereof. 
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I. INTRODUCTION 

During the last two decades the Hubbard Hamiltonian 
has become the model for strongly correlated particles, 
both fermions and bosons, and nowadays its investigation 
does not need to be justified. A great variety of phenom- 
ena have been interpreted in terms of the Hubbard model, 
such as ferromagnetism, antiferromagnetism, bond alter- 
nation, spin liquids, superconductivity and even charge 
stripes. Due to the complexity of the theoretical prob- 
lem the analysis is typically based either on approxi- 
mate treatments - mean-field theory, perturbative ex- 
pansions, variational wave functions - or on "numerically 
exact" methods - diagonalization of the Hamiltonian for 
small system sizes, quantum Monte Carlo, Density Ma- 
trix Renormalization Group. Each method has its own 
advantages, but also its drawbacks. Clearly one has to 
examine the various results with a critical mind because 
often the approximations are too crude and may thus 
overestimate certain ordering tendencies. A famous ex- 
ample is the Stoner theory of ferromagnetism. If applied 
to the simple Hubbard model on the square lattice at 
half filling, this theory predicts a ferromagnetic ground 
state for any (positive) value of the Hubbard parame- 
ter U. In contrast, detailed variational calculations show 
that there is no ferromagnetism at half filling (where the 
ground state is antiferromagnetic) and close to half filling 
ferromagnetic ordering could at most occur for extremely 
large values of U. 

While ferromagnetic long-range order fades away if 
methods are used that are more accurate than Stoner 
theory, the situation is radically different in the context 
of superconductivity. In fact, within BCS theory the 
repulsive Hubbard model is not superconducting. The 
on-site repulsion clearly is detrimental to on-site pair- 
ing, which can be discarded right away. For other types 
of pairing, such as extended s-wave, p-wave or d-wave, 
the on-site repulsion gives zero contribution to the con- 
densation energy in BCS mean-field theory. Thus, if the 
Hubbard model has a superconducting ground state, then 
the order parameter has to be stabilized by electron cor- 
relations. A pictorial view of such an unusual supercon- 
ductor has been given by Anderson in his theory of res- 
onating valence bonds where singlet bonds of the parent 
antiferromagnetic insulator turn into charged supercon- 
ducting pairs upon doping 1 . Detailed recent studies for 



small, intermediate and large values of U confirm that 
the repulsive Hubbard model on a square lattice does ex- 
hibit superconductivity with <i-wave symmetry for elec- 
tron densities close to (but not equal to) 1. 



II. FROM SMALL TO LARGE U 

We consider the Hubbard Hamiltonian 
H = Ho + UD , 

where 

H = - ^2u 3 c\ a c J<7 



(1) 



(2) 



describes hopping over the sites of a square lattice, with 
the usual fermion creation and annihilation operators, cl 
and Ciu, respectively, and 



with ru 



(3) 



measures the number of doubly occupied sites. Only hop- 
ping between nearest (ty = i) and next-nearest neigh- 
bours (Uj = t') will be considered. The bare spectrum, 

£k = — 2f (cos k x + cos k y ) — At' cos k x cos k y , (4) 

has a bandwidth W — 8t. We discuss the three regimes 
of small U (U < W), intermediate U (U « W) and 
large U (U ^> W) separately because they are usually 
approached with rather different techniques. 



A. Small U 

The question of Cooper pairing in a fermionic system 
with purely repulsive interactions has been addressed a 
long time ago by Kohn and Luttinger— using a pertur- 
bative calculation of the effective two-particle vertex in 
powers of the coupling strength. Their estimate for the 
critical temperature in the case of 3 He was deceptively 
small, but it was argued that in electronic systems with 
an appropriate band structure the situation may be more 
favorable^. In the case of the Hubbard model on the 
square lattice close to half filling, divergences preclude 
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the application of naive perturbation theory, but the de- 
velopment of the renormalization-group approach for in- 
teracting fermions during the early nineties paved the 
way for studying the effective two-particle vertex also 
in this case. Several groups used this technique, which 
consists in tracing out high-energy states in incremen- 
tal steps and at the same time incorporating their ef- 
fects within a renormalized vertex, which is a function of 
momenta^—. These studies differ in technical details, but 
not in the main message, namely that at and very close to 
half filling the ground state of the two-dimensional Hub- 
bard model is antiferromagnetic while for higher but not 
too high doping d-wave superconductivity prevails. It is 
not easy to provide accurate quantitative predictions for 
gap sizes, critical temperatures or order parameters us- 
ing this approach, but for U of the order of i/2 these 
quantities are not ridiculously small. 

B. The large U limit 

In the large U limit, where double occupancy is sup- 
pressed, the Hubbard model can be replaced by the t-J 
model close to half filling 7 - , defined by the Hamiltonian 

-fft-J = - ^J C L C JCT + ^ JijSi ■ Sj , (5) 

ija ij 

where Si are spin i operators and J t j = 2t?j/U. This 
Hamiltonian is supposed to act in the reduced space of 
configurations without doubly occupied sites. It is im- 
portant to mention that in general additional three-site 
terms exist which are expected to be negligible only very 
close to half filling. Further away from half filling these 
terms should be kept if one wants to study the large U 
limit of the Hubbard model. A great variety of techniques 
have been applied to the t-J model^^. Here we just men- 
tion a simple variational ansatz for the ground state, the 
"resonating valence bond state" (RVB) of Anderson^, 

l*RVB) =II( 1 - n W)l*o}, (6) 

i 

where \^f ) represents the filled Fermi sea or a BCS sin- 
glet superconductor. Variational Monte Carlo results for 
d-wave pairing with a gap function 

A(k) = A (cosfc :E — cos k y ) (7) 

yield a large region of doping, up to about 40%, where 
this state is preferred^. The gap parameter Ao is found 
to be largest at half filling, of the order of t, at the same 
time the pair amplitude goes to zero when approaching 
half filling. This unusual behavior has been associated 
with the pseudogap of underdoped cuprates^. 

C. Cuprates, a case for intermediate U 

The application of the Hubbard model to layered 
cuprates is a priori not obvious, but several experiments 



indicate that this is a good starting point. Photoemission 
data exhibit a single band crossing the Fermi surface, at 
least in optimally and overdoped samples, and the shape 
of the Fermi surface can be fitted rather welli^ using the 
tight-binding band structure (|4]) with t' « — 0.3i. Inelas- 
tic neutron scattering for the parent antiferromagnetic 
insulator can be interpreted in terms of the Heisenberg 
model, provided that ring exchange is included^. This 
term appears naturally in the Hubbard model if the large 
U expansion is carried out to higher orders than in Eq. 
([5]) . The comparison between theoretical and experimen- 
tal dispersion curves for magnetic excitations givesi^ an 
estimate of U ~W. The analysis of the neutron data is 
consistent with photemission experiments on the antifer- 
romagnetic insulator, where the momentum distribution 
n(k) shows a pronounced k-dependencei^, in contrast to 
the t-J model where n(k) = 1. We conclude that U 
is of the order of the bandwidth W where neither the 
functional renormalization group, valid for U <C W, nor 
the t-J model, representing the U 3> W regime, can be 
trusted. 



D. Variational results for intermediate U 

Numerical approaches, such as quantum Monte 
CarloiS or quantum cluster method a 17 ' 18 are widely used 
to unravel the secrets of the Hubbard model, especially 
for intermediate values of U. These techniques have cer- 
tain advantages, for instance to be rather unbiased, but 
they have also their problems, such as the trouble with 
the fermionic sign and/or the restriction to small cluster 
sizes. An alternative route is the use of variational wave 
functions. The RVB state (|6|) is clearly not appropriate 
for U of the order of the bandwidth, where doubly oc- 
cupied sites cannot be excluded. Many other proposals 
have appeared in the literature, but we restrict ourselves 
to one of them, defined as 

|*GB)=e- hfl °/ t e-^|* >. (8) 

The operator e _gD partially suppresses double occupancy 
for g > 0, while c^ 1111 "/* promotes both hole motion and 
kinetic exchange. The limit h — > leads to the Gutzwillcr 
ansatz. In the large U limit and close to half filling the 
variational ansatz (|5]) for the Hubbard model is equiv- 
alent to the RVB state (U5J) applied to the t-J models, 
but we should keep in mind that the state ([8]) is a rather 
poor approximation for U 3> W 19 i 20 . 

The results obtained with the variational ground state 
(|8l) for U = 8t, t' = — 0.3t can be summarized as 
follows 21 . Superconductivity occurs in a limited doping 
range, namely for electron densities 0.75 < n < 1 on the 
hole-doped side and 1.05 < n < 1.2 for electron dop- 
ing. The two endpoints 0.75 and 1.2 correspond to the 
two specific densities beyond which the antiferromagnetic 
zone boundary does no longer cut the Fermi surface (the 
"hot spots" have disappeared). This strongly hints at 
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pairing induced by antiferromagnetic fluctuations. Very 
close to half filling the magnetic correlations become so 
strong that they suppress both the superconducting gap 
parameter and the Cooper pair amplitude. There are 
qualitative differences between hole and electron doping. 
In the former case the energy gain comes clearly from 
the kinetic energy, while in the latter case the more con- 
ventional gain in potential energy is found. Moreover 
the condensation energy for the hole-doped side is much 
larger than that for electron doping. 

The variational results obtained with the ansatz ((SJ) arc 
in good agreement with quantum cluster calculations 22 
and experimental phase diagrams 2 ^. A quantum Monte 
Carlo study by Aimi and Imada 24 , often quoted as 
counter-argument against superconductivity in the repul- 
sive Hubbard model, is also consistent with these varia- 
tional results, not however with the predictions of the 
RVB trial state ©. 



III. CONCLUDING REMARKS 



An important issue in the discussions of layered 
cuprates is the nature of the "pseudogap phase". In 
contrast to the RVB theor y 11 ! 12 the variational results 
discussed above did not show an energy scale steadily in- 
creasing when approaching half filling. This can be due 
to the fact that competing instabilities, such as charge 
stripes, incommensurate spin-density waves or circulat- 
ing currents have not been included in the ansatz. In any 
case, the superconducting gap, which was found to de- 
crease when approaching half filling, is more naturally as- 
sociated with an onset temperature for superconducting 
fluctuations than with the temperature T* below which 
the pseudogap appears. 
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